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ABSTRACT 

Searches for companions of brown dwarfs by direct imaging mainly probe orbital separations greater than 3-10 AU. On the other 
hand, previous radial velocity surveys of brown dwarfs are mainly sensitive to separations smaller than 0.6 AU. It has been speculated 
that the peak of the separation distribution of brown dwarf binaries lies right in the unprobed range. This work extends high-precision 
radial velocity surveys of brown dwarfs for the first time out to 3 AU. Based on more than six years UVES/VUI spectroscopy the 
binary frequency of brown dwarfs and (very) low-mass stars (M4.25-M8) in ChamaeleonI was determined: 18^[2 % for the whole 
sample and 10;^^** % for the subsample of ten brown dwarfs and very low-mass stars (M< 0.1 Mq). Two spectroscopic binaries were 
confirmed, the brown dwarf candidate ChaHo- 8 (previously discovered by Joergens & Miiller) and the low-mass star CHXR74. Since 
their orbital separations appear to be 1 AU or greater, the binary frequency at <1 AU might be less than 10%. Now for the first time 
companion searches of (young) brown dwarfs cover the whole orbital separation range, and the following observational constraints for 
models of brown dwarf formation can be derived: (i) the frequency of brown dwarf and very low-mass stellar binaries at <3 AU does 
not significantly exceed that at >3 AU; i.e. direct imaging surveys do not miss a significant fraction of brown dwarf binaries; (ii) the 
overall binary frequency of brown dwarfs and very low-mass stars is 10-30%; (iii) the decline in the separation distribution of brown 
dwarfs towards smaller separations seems to occur between 1 and 3 AU; (iv) the observed continuous decrease in the binary frequency 
from the stellar to the substellar regime is confirmed at <3 AU providing further evidence of a continuous formation mechanism from 
low-mass stars to brown dwarfs. 

Key words. binaries: spectroscopic — planetary systems — stars: individual ([NC98] Cha HA 1, [NC98] Cha HA 2, 
[NC98] Cha HA 3,[NC98] Cha HA 4, [NC98] Cha HA 5, [NC98] Cha HA 6, [NC98] Cha HA 7, [NC98] Cha HA 8, 
[NC98] Cha HA 12, CHXR74, CHXR76 (B34)) — stars: low-mass, brown dwarfs — stars: pre-main sequence — techniques: radial 
velocities 



1. Introduction 

I Searching for companions to brown dwarfs (BDs) is essential 
for understanding BD formation, for which no widely accepted 
' model exits (e.g. recent review by Luhman et al. 2007). A priori, 
I one might expect that BDs form in the same manner as stars by 
■ the direct collapse and fragmentation of molecular clouds, only 
on a much smaller scale. Standard cloud fragmentation seems 
to have difficulty in making BDs (Boss 2001; Bate, Bonnell & 
Bromm 2003). One possible explanation is that the simulations 
lack an important piece of physics; e.g., supersonic magneto- 
turbulence can create local high densities and facilitate the for- 
mation of BDs (Padoan & Nordlund 2004). Alternatively, BDs 
are formed when cloud fragmentation is modified by an addi- 
tional process that prematurely halts accretion during the proto- 
stellar stage, such as dynamical ejection out of the dense gaseous 
environment (e.g. Reipurth & Clarke 2001; Boss 2001; Bate, 
Bonnell & Bromm 2002; Sterzik & Durisen 2003; Umbreit et 
al. 2005) or photoevaporation by ionizing radiation from nearby 
massive stars (Kroupa & Bouvier 2003; Whitworth & Zinnecker 
2004). A fourth, recently studied scenario foresees BD formation 
through instabilities in the outer (>100AU) disk around a star 



* Based on observations collected at the European Southern 
Observatory, Chile in program 75.C-0851(C), 76.C-0847(A), 77.C- 
0831(A+D). 



and release into the field by e.g. secular perturbation (Goodwin 
& Whitworth 2007; Stamatellos et al. 2007). 

The frequency and properties of BDs in multiple systems are 
fundamental parameters in these models, e.g. embryo-ejection 
scenarios for BD formation predict only a few BD binaries in 
preferentially close orbits, while isolated fragmentation scenar- 
ios are expected to generate BD binaries with similar (scaled- 
down) properties, as stellar binaries have. However, the binary 
properties of BDs are poorly constrained for close separations. 
Most of the current surveys for companions to BDs and very 
low-mass stars (VLMS, M < 0.1 Mq) are done by direct (adap- 
tive optics or HST) imaging (e.g. Reid et al. 2001; Bouy et 
al. 2003; Burgasser et al. 2003; Gizis et al. 2003; Close et al. 
2003; see also Burgasser et al. 2007) and are not sensitive to 
either close separations (a < 2-3 AU and a < 10-15 AU for 
the field and clusters, respectively) nor to low-mass ratio sys- 
tems (qsM2/Mi <0.6). These surveys find a lower frequency 
(10-30%) of BD/VLM binaries compared to stars, a separation 
distribution with a peak around 3-10 AU, and a preference for 
equal-mass components with 77% having a mass ratio q>0.8. 
The observed peak of the separation distribution is close to the 
incompleteness limit; thus, we do not know if there is a signif- 
icant number of still undetected, very close (<3 AU) BD/VLM 
binaries. It has been suggested that BD/VLM binaries with a 
separation <3 AU are as frequent or even more frequent than 
those at >3AU (Pinfield et al 2003; Maxted & Jeffries 2005; 
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Burgasser et al. 2007). Hence, the peak of the BD/VLM binary 
separation distribution could he below 3 AU. Furthermore, while 
it has been argued that the higher frequency of systems with 
mass ratios close to unity appears to be a real feature of very 
low-mass binaries, it would be, nevertheless, desirable to actu- 
ally probe mass ratios lower than 0.6. Given these limits to the 
current observational data, the question arises as to whether our 
current picture of BD/VLM binaries is complete, or, whether we 
miss a significant fraction of very close and/or low-mass ratio 
systems. 

Spectroscopic monitoring of radial velocity (RV) variations 
provides a means to detect very close binaries. Given sufficient 
RV precision, the detection of low-mass ratio systems is fea- 
sible down to planetary masses (cf. Joergens 2006a; Joergens 
& Miiller 2007). While several spectroscopic surveys for com- 
panions of BD/VLMS in young clusters (Joergens & Guenther 
2001; Kenyon et al. 2005; Joergens 2006a; Kurosawa, Harries 
& Littlefair 2006; Prato 2007a; Maxted et al. 2008) and in 
the field (Reid et al. 2002; Guenther & Wuchterl 2003; Basri 
& Reiners 2006) have been started in recent years, data sam- 
pling is sparse in most cases and the number of confirmed 
close companions to BD/VLMS is still small. To date, there 
are four spectroscopic BD/VLM binaries known for which an 
orbital solution has been determined: PP115 (Basri & Martin 
1999), GJ569Bab (Zapatero Osorio et al 2004; Simon, Bender 
& Prato 2006), 2MASS J05352184-0546085 (Stassun, Mathieu 
& Valenti 2006), and ChaHaS (Joergens & Muller 2007). One 
of them, ChaHaS, was discovered (Joergens & Muller 2007) 
within an RV survey for (planetary and BD) companions of 
very young BD/VLMS in the Chamaeleon I star-forming region, 
which is being caiTied out with UVES at the VLT (Joergens 
& Guenther 2001; Joergens 2006a). ChaHaS is presumably 
the lowest mass RV companion detected so far in a close orbit 
around a BD/VLMS and only the second very young BD/VLM 
spectroscopic binary. That survey in Cha I provided the first de- 
termination of the binary frequency of very young BD/VLMS 
at close separations and the indication that the fraction of un- 
detected BD binaries at close separations must be small. These 
results were based on probing orbital separations below 0.3 AU. 
The present paper reports follow-up spectroscopic observations 
and RV measurements of the targets of that survey. This signif- 
icantly enhances the separation range probed. For the first time 
the binary frequency of (very young) BD/VLMS at separations 
of 3 AU and smaller is determined. 

The paper is organized as follows. Section|2] describes the 
spectroscopic observations and RV measurements, in Sect.[3]the 
identification of RV variable objects is presented, in Sect. |4] the 
observed binary fraction is derived and the probed separation 
range assessed, Sect.|5]gives details for the spectroscopic system 
CHXR74, and finally, Sect.|6]and Sect.|7] provide a discussion 
and conclusions. 



2. Observations and radial velocities 

Spectroscopic follow-up observations of BDs and (very) low- 
mass stars (VLMSs) in Cha I were carried out between October 
2005 and January 2006 with the Uv- Visual Echelle Spectrograph 
(UVES) attached to the VLT S.2m KUEYEN telescope at a 
spectral resolution A/AA of 40000 in the red optical wave- 
length regime. In program 76.C-0S47(A), UVES spectra were 
taken of all objects that were previously observed by Joergens 
(2006a) and that still lacked re-observation after more than 
one year (Cha Hal, Cha Ha 2, Cha Ha 3, Cha Ha 5, Cha Ha 6, 



ChaHalfl and B3^. Though Cha Ha 7 was also a tar- 
get of this run, observations failed for it because of wrong 
pointing of the telescope in service mode. In programs 75. C- 
0S51(C) and 77.C-0S31(A-hD), follow-up spectra of the low- 
mass star CHXR 74 were taken, which was previously revealed 
as an RV variable source (Joergens 2006a). Determinations of 
basic (sub)stellar parameters of the objects can be found in 
Comeron, Rieke & Neuhauser (1999), Comeron, Neuhauser 
& Kaas (2000), and Luhman (2004, 2007). It is noted that, 
while the cuiTent knowledge of the Cha I cloud is complete to 
0.01 Mq (Luhman 2007), this sample is biased towards brighter 
BD/VLMS. The reason for this is first, that the sample was cho- 
sen from late-M type objects that were known at the time of the 
start of this survey in 2000 (Comeron et al. 2000) from shallower 
surveys, and, second that three of the faintest objects known at 
that time, Cha Ha 9, Cha Ha 10, and Cha Ha 11, were not in- 
cluded. 

After standard CCD reduction of the data, one-dimensional 
spectra were extracted and their wavelength scale established 
in a first step by means of ThAr lamp exposures. The RVs 
were measured from these spectra based on a cross-coiTelation 
technique employing telluric lines for the wavelength calibra- 
tion. Details on the data reduction and analysis can be found in 
Joergens (2006a). Each RV determination is based on two con- 
secutive single spectra that provide two independent measure- 
ments of the RV (an exception here are the 2006 observations of 
CHXR 74). This allows an estimation of the eiTor of the relative 
RVs based on the sample standard deviation for two such data 
points. The new RV measurements are summarized in Table[T] 
These estimated errors of the relative RVs range between 30 
and 100ms ' (Table[T]). The absolute RVs are subject to an ad- 
ditional error of about 400 ms"' because of the uncertainty of 
the velocity zero point (cf. Joergens 2006a). Figures[T]-@]display 
the results. The studied objects have a mean RV of 16.0 km s"' 
in agreement with membership in the Cha I star-forming cloud 
(the surrounding molecular gas has a velocity of 15.4kms"', 
Mizuno et al. 1999). Their RV dispersion is relatively small: 
the RV dispersion measured in terms of standard deviation is 
O.Skms"' and the RVs cover a total range of 2.2 km s"', con- 
sistent with previous kinematic measurements of the objects 
(Joergens 2006b). In particular, no outliers are found. 



3. RV Variability 

There are different approaches to identifying significant vari- 
ability in datasets of repeated RV measurements. Prato (2007b) 
classifies an object as RV variable if the root-mean-square (rms) 
scatter of its repeated RV measurements are significantly (five 
times) above the average rms scatter of the whole sample, i.e. as- 
suming that most objects are not intrinsically RV variable. This 
method registers external noise, like RV noise caused by activ- 
ity, by calculating the average rms scatter of the sample. Other 
authors compute the probability with which observed RV values 
of an object with individual eiTor estimates can be described by 
a constant function in a test (e.g. Maxted & Jeffries 2005). 
Objects for which the probability for being RV constant is less 
than a chosen threshold are classified as RV variable sources. 
The identification of RV variability and, therefore, of spectro- 
scopic binaries by this method depends crucially on the individ- 
ual measurement eiTors. 



' Simbad names: [NC98] Cha HA 1, [NC98] Cha HA 2, etc. 
- Simbad name: CHXR 76 
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For the identification of RV variable objects among 
BD/(V)LMS in Cha I, the second method is appHed and the 
probabiHty p for each object being RV constant is computed. 
The constant function, which is fitted in this test, is given by 
the weighted mean of the individual RV values. The threshold 
of tolerance for error in the;\f^ test is set to 0.1% (p < 10^^), 
i.e. a possible detection of RV variability has a significance of 
at least 3.3 cr. The evaluated data set for BDs and (very) low- 
mass stars in Cha I includes the new RV measurements presented 
here (Table[T]i as well as previous ones by Joergens (2006a) and 
Joergens & Miiller (2007). Table|2]summarizes the relevant char- 
acteristics of these data for each studied object (half peak-to- 
peak RV differences, the weighted mean RV, the error of this 
mean, which takes into account an uncertainty of 400 ms"' for 
the zero point of the velocity, the number of observations, the 
minimum time separation between two RV measurements, and 
the covered time base). For the variability classification, in ad- 
dition to the listed internal RV errors (in most cases standard 
deviation of two consecutive measurements) a minimum error 
of 0.3 km s ' has been applied to account for fluctuations in the 
standard deviation for a small number of measurements and/or 
for systematic errors. This minimum error is derived by an es- 
timate of the systematic errors based on the night-to-night rms 
scatter of all BD/VLMS in the sample of 0.32 km s ' on aver- 
age (i.e. without making any assumptions on the internal errors 
or the variability status), which is in good agreement with their 
average internal error of 0.29kms The calculated proba- 
biUties are listed in TablelH Both Cha Ha 8 and CHXR74 have 
values below the chosen threshold {p < 10"^), hence, are iden- 
tified as RV variables. For both, Cha He 8, which was recently 
shown to have a BD companion (Joergens 2006a; Joergens & 
Muller 2007) and for CHXR74, for which a longterm RV trend 
was observed (Joergens 2006a), these are confirmations of previ- 
ous findings of spectroscopic companions. For CHXR74, there 
are two epochs of RV data added in this work providing further 
constraints on the orbit of the companion. See Sect.|5]for more 
details on CHXR74. 

Figure|5] shows the distribution of probabilities. For the 
RV constant objects, it resembles the expected distribution of 
- log p for random fluctuations alone normalized to the num- 
ber of constant objects (dashed line). This shows that the ap- 
plied error estimates are reliable and, in particular, that they have 
not been underestimated. The RV variable objects, Cha Ha 8 and 
CHXR 74, appear in the rightmost bin of the histogram in Fig.|5] 
They are clearly separated from the nonvariable group. 

The performed variability classification by ax^ test is based 
on the assumption that the errors are independent of the data. 
Applying a minimum error given by the average scatter of the 
measurements, as done here, accounts for non-Gaussian errors 
and partly ensures such an independence. To further increase the 
confidence, the result has been double checked by means of an 
analysis of variance (ANOVA) test, which is a generalization of 
the student's t-test to more than two measurements and which is 
independent of error estimates. The F probabilities were calcu- 
lated using the individual RV measurements, i.e. for RV values 
that are based on two consecutive RV measurements, the two 
individual RV values were used. It is found that the outcome of 
this analysis of variance is basically consistent with that of the;\f^ 
test; in particular, the same objects were found as significantly 
variable sources. 

One can use the presented data to assess the level of RV noise 
for very young BD/VLMS (in Cha I) by quantifying the night-to- 
night rms scatter of only the constant BD/VLMS in the sample, 
which is 0.31 kms"' on average. 



Another way to look at the results of this survey is to con- 
sider the recorded RV differences. Figure[T3]displays half peak- 
to-peak RV differences ARV versus the object mass, as adopted 
from Comeron et al. (1999, 2000fl Objects classified as RV 
constant have ARV <1 kms"', and the majority of them have 
even smaller ARV below 0.3 kms '. The two RV variables, on 
the other hand, exhibit RV differences ARV above 1.4kms 
Though clearly separated from the vast majority of the RV con- 
stant objects in this diagram, it is remarkable that RV variability 
occurs with such relatively small amplitudes. This is discussed 
further in Sect.|6] 

It is noted that, because of the small variability amplitudes, 
none of the objects would have been identified as RV variable 
source when applying the criterion of Prato (2007b), namely 
to require that the rms scatter has to be five times greater than 
the average rms scatter of the whole sample to classify an ob- 
ject as RV variable: The rms scatter of the RV variable ob- 
jects, Cha Ha 8 and CHXR 74, is 1.0 and 1.5 kms respec- 
tively, while the average rms scatter of the whole sample is 
0.6 km s"' and that of the subsample of BD/VLMS is 0.5 km s"' . 
It is further noted that, if the intrinsic errors of the RV mea- 
surements would have been as large as in the survey of Basri & 
Reiners (2006, 1.4 km s"' on average), also none of the objects 
would have been identified as RV variable. 

In a similar survey but with more restricted separation 
range, Kurosawa, Harries & Littlefair (2006) find four spec- 
troscopic BD/VLM binary candidates in USco and p Oph. The 
recorded half peak-to-peak RV differences for them (ARV = 0.2- 
0.5 km s"') are significantly less than for the variable BD/VLMS 
in Chal (l-2kms"'). Figure[T4lshows ARV plotted versus mass 
for the data from Kurosawa et al. (2006). It can be seen here 
that the four RV variable objects (marked in the figure with their 
names) cannot be distinguished from the RV constant objects 
based on larger RV differences. (The average intrinsic error of 
the whole sample is 0.42 kms"', five objects in the sample have 
RV errors <0.15kms"', four of them are found to be signifi- 
cantly variable.) While it has to be further investigated whether it 
can be excluded that the detected variability is caused by system- 
atic errors, if we interpret them as RV semi-amplitudes caused 
by companions, these companions would have masses M2 sin / in 
the planetary regime (1-2 Jupiter masses) assuming a separation 
of <0.1AU. 



4. Binary fraction and probed orbital separations 

Based on the statistical analysis in the previous section, two 
objects (Cha Ha 8, CHXR74) among a sample of eleven BDs 
and (very) low-mass stars (spectral type M4.25 and later; M< 
0.25 Mq) are found with significant RV variations and are, there- 
fore, considered as spectroscopic binary systems. The uncer- 
tainty of the measured binary fraction is estimated by search- 
ing for the fraction of binaries at which the binomial distribution 
PB(x,n,p) for X positive events in n trials falls off to 1/e of its 
maximum value (cf. Basri & Reiners 2006). The observed bi- 
nary fraction of the whole sample is 18;^j2 % (2/11). When con- 
sidering only the BD/VLMS in the sample, i.e. the ten objects 
with M< 0.1 Mq, we observe one RV variable (Cha Ha 8), i.e. a 
BD/VLM binary fraction of lO+i*^ % (1/10). 



^ Luhman (2004, 2007) derives systematically higher effective tem- 
peratures and luminosities for the objects leading to higher mass 
estimates when compared with evolutionary tracks, e.g. 0.1 for 
Cha Ha 8 instead of 0.07 Mq. 
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The orbital separation range that can be probed by an RV 
survey for companions depends mainly on the timely sampling 
of the measurements and on their uncertainties. To assess the 
detection efficiency as a function of the orbital separation of this 
RV survey, a Monte-Carlo simulation of RV measurements for 
randomly selected binaries from a model was performed using 
the time spacing and error estimates as given by the real obser- 
vations and again applying a minimum error of 0.3 km s"^ For 
each simulated binary, it is then checked whether it would have 
been detected as RV variable in this survey, i.e. if its probabil- 
ity for being RV constant is lower than the chosen threshold of 
p < 10"-' (cf. Sect.|3]l. A similar but more sophisticated simu- 
lation has been performed by Maxted & Jeffries (2005) and ap- 
plied in a simplified version by other authors (e.g. Kurosawa et 
al. 2006). The procedure is described in these works and the fol- 
lowing summarizes only the basic steps and assumptions applied 
here. For each object of the sample, the relative RV errors, num- 
ber of RV data points, and their temporal separation as given by 
the real observations are considered. Furthermore, the primary 
mass estimates from Comeron et al. (1999, 2000) are used. For 
each semi-major axis in the range between 0.006 and 10 AU, 10^' 
binaries are simulated by (i) assuming a flat mass ratio distribu- 
tion between q = Mil Mi =0.2 and 1.0, i.e. selecting the mass 
ratio randomly from this range; (ii) assuming circular orbits, i.e. 
setting the eccentricity to zero and not considering the longitude 
of periastron; (iii) selecting the orbital inclination / randomly 
from a uniform cos / distribution; and by (iv) selecting the or- 
bital phase for the first simulated RV randomly between and 
1 . These assumptions are justified by the results of Maxted & 
Jeffries (2005), who find that the detection probability is mainly 
insensitive to the distribution of mass ratio and eccentricity. The 
minimum considered separation in this simulation is 0.006 AU. 
This allows a distance between central object and companion of 
more than 2.4 times the radius (Roche criterion for stable orbits) 
and assumes a size of half a solar radius. (Very young BDs have 
relatively large radii, e.g. the radius of a BD in Orion with mass 
0.034 Mq was measured to 0.51 7?o=0.0024 AU, Stassun et al. 
2006.) For each binary selected in the described way, RV mea- 
surements are simulated by calculating the RV for the first or- 
bital phase, evolving the orbit based on the given time sequence 
of the real observations, and calculating the subsequent RV mea- 
surements. These simulated RV data, along with the RV errors 
of the real observations, are then subject to the same RV vari- 
ability check as performed for the observations. Finally, the de- 
tection probability is given by the number of detected binaries 
compared to the total number of 10^ simulated binaries for each 
semi-major axis. It is noted that line-blending is not taken into 
account in this simulation. This effect is most relevant for nearly 
equal-mass systems. Maxted et al. (2008) show that neglecting 
line-blending can lead to an overestimate of the detection effi- 
ciency of up to 15%. This has to be kept in mind when interpret- 
ing the detection probabilities calculated here and in the frame- 
work of other RV surveys of BD/VLMS that neglect this effect 
(e.g. Kurosawa et al. 2006; Basri & Reiners 2006). 

Figure|6] shows the results of the Monte-Carlo simulation 
for all objects in the order of increasing sensitivity for larger 
separations. The majority of objects in the sample (ChaHa 1, 
ChaHa2, ChaHa 3, ChaHa 5, ChaHa 6, ChaHa 12) was ob- 
served with a temporal spacing of about 20 days and re- 
observation after more than about 2100 days (5.8 years). Four 
objects (ChaHa 4, ChaHa 8, B 34, CHXR74) were observed 
with a higher cadence. For them, between 6 and 13 indepen- 
dent RV data points were obtained covering a total time base 
between 700 and 2540 days (2-7 years). As can be seen from 



Fig.|6] the range of orbital separations that can be probed with a 
50% or higher detection probability is <1 AU for ChaHa 1 and 
Cha Ha 3, < 3 AU for Cha Ha 2, Cha Ha 4, Cha Ha 5, Cha Ha 6, 
ChaHa 8, and ChaHa 12, <4AU for B34, and <6AU for 
CHXR74, respectively. For Cha Ha 7 (Fig.|51 top panel), the 
probed separation range is significantly narrower (<0. 2 AU) than 
for the other objects because scheduled re-observations failed 
(cf. Sect.|2]i. For this BD, the total time base of its RV mea- 
surements is 20 days. On average, the semi-major axis range of 
<3 AU was probed with a 50% or higher detection rate. 

The detection probability can be translated into a rate of bi- 
nary systems present in the sample but missed in the survey. 
However, the absolute number of missed binaries in this survey 
cannot be quantified because the orbital period distribution is un- 
known for BD/VLMS. Therefore, no corrected binary fraction is 
derived based on the observed binary fraction. 

The maximum separation to which an RV survey is sensitive 
depends not only on the total covered time but, in addition, also 
crucially on the measurement errors. For example, Cha Ha 5 has 
been observed with almost exactly the same time pattern as e.g. 
Cha Hal and Cha Ha 3; however, the separation range probed 
with 50% detection efficiency is <3 AU for Cha Ha 5 compared 
to < 1 AU for Cha Ha 1 and Cha Ha 3, respectively. The only dif- 
ference here is the RV error estimates, which are on average 
0.4kms ' forChaHa5 and 0.5kms"' for the other two objects 
(applying a minimum error of 0.3 km s"'), respectively. 

Figure|2] shows that assuming a fixed value for the orbital 
inclination in the Monte-Carlo simulation, as applied by Basri 
& Reiners (2006), leads to a slight overestimation of the detec- 
tion efficiency for larger orbital separations compared to the re- 
alistic case of random orientation. In Fig.|8] it is demonstrated 
that third-epoch observations between the observations with the 
largest time difference can have a significant effect on the proba- 
bility distribution. While such intermediate third-epoch observa- 
tions do not influence the maximum probed separation, they can 
increase the detection rate at smaller separations, here by up to 
20%. 

In the following, the detection efficiency vs. orbital separa- 
tion of several current RV surveys of BD/VLMS (Guenther & 
Wuchterl 2003; Kurosawa et al. 2006; Basri & Reiners 2006; 
this work) is compared in a uniform way. For this purpose, the 
detection probability is calculated (i) by taking into account the 
average measurement uncertainties and time sequences of the 
individual surveys; (ii) by assuming circular orbits, random ori- 
entation, random starting orbital phase, a mass ratio chosen ran- 
domly between 0.2 and 1 ; and (iii) by applying a threshold for 
variability detection of p < 10 Figures|9HT0l show the result. 
The orbital separation range that can be probed with 50% detec- 
tion efficiency is <3 AU (subsample of BD/VLMS in this work), 
<0.6 AU (Guenther & Wuchterl 2003), <0.3 AU (Kurosawa et al. 
2006), and <0.2AU (Basri & Reiners 2006), respectively. The 
difference in the maximum orbital separations between the sur- 
veys by Kurosawa et al. (2006), Guenther & Wuchterl (2003), 
and this work can be largely attributed to differences in the cov- 
ered time base. While the distribution of the detection efficiency 
in the survey by Basri & Reiners (2006) is non-zero for separa- 
tions greater than 1 AU due to the relatively large time base, 
it is decreasing more rapidly compared to the other surveys. 
This is mainly due to the moderate RV precision of this sur- 
vey, which makes it less sensitive to low-mass ratio systems. 
The differences of this calculated detection probability curve for 
the survey by Basri & Reiners (2006) compared to the one pre- 
sented by the authors themselves, who find a sensitivity up to 
6 AU, are the assumptions concerning inclination (random ori- 
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entation instead of fixed inclination), mass ratio {q selected ran- 
domly from a uniform distribution between 0.2 and 1 instead of 
^=0.5), and threshold for variability (p < 10"^, i.e. 3.3cr instead 
of p < lO"''^"*^, i.e. 2cr). This, in addition to neglecting differ- 
ences in the RV errors of different surveys, leads to a deviating 
picture in Fig. 4 of Basri & Reiners (2006). 

The BD/VLM binaries with a mass ratio close to unity cause 
a stronger RV signal compared to low-mass ratio systems and 
are, therefore, easier detected at larger separations. To quantify 
this, the sensitivity of the considered surveys has been calcu- 
lated for binaries with mass ratios between 0.8 and 1 . In this case 
iq >0.8), a 50% detection efficiency can be read off from Figs.fTTI 
and fT2l for the following separation ranges: <5 AU (subsample 
of BD/VLMS in this work), <2 AU (Basri & Reiners 2006), 
<0.8 AU (Guenther & Wuchterl 2003), and <0.4 AU (Kurosawa 
et al. 2006). In this mass-ratio regime, neglecting line-blending 
in the performed simulation is relevant (see above). At separa- 
tions where the effects of line blending in the q range 0.8-1.0 are 
significant, these detection probabilities may be significantly re- 
duced. A comprehensive analysis of these effects is beyond the 
scope of this work. 

For eleven objects in the survey by Basri & Reiners (2006), 
there are a third or more epoch observations with a temporal sep- 
aration of more than one day. Such intermediate observations 
can positively influence the detection efficiency, as shown be- 
fore. However, a simulation assuming an average observing pat- 
tern of fofej=[0,500d,1650d] for this survey does not change the 
calculated detection probability. 

5. Low-mass, long-period spectroscopic binary 
CHXR74 

Significant RV variability indicating the presence of a spec- 
troscopic companion orbiting the very young low-mass star 
CHXR74 (M4.25-M4.5,~O.2-O.25M0, Comerdn et al. 2000; 
Luhman 2007) has been akeady found by Joergens (2006a). 
The new RV measurements of CHXR74 in 2005 and 2006 
(Table[T]) support this finding. The RV of this star is, apart from 
some short-term scatter, increasing almost monotonically be- 
tween 2000 and 2006 (Fig.lH bottom panel) with a recorded 
peak-to-peak difference of 4.8 km s '. The data indicate a rel- 
atively long orbital period of at least twelve years, hence, an or- 
bital separation greater than about 3 AU. This can be used to 
derive a very rough estimate of the minimum companion mass 
M2 sin / of ~80Mji,p. Given the relatively long orbital period for 
a spectroscopic system and its youth, CHXR 74 might be one 
of the first very young spectroscopic binaries that can be spa- 
tially resolved. A binary with a separation of more than 3 AU 
(20 milli arc sec at the distance of the target) might be spa- 
tially distinguished from the primary by high-resolution imag- 
ing, e.g. with NACO/VLT. If so, CHXR 74 might provide dy- 
namical mass determinations of its components and, therefore, 
valuable constraints for the early evolution at the low-mass end 
of the stellar mass spectrum. 

6. Summary and discussion 

This paper presents an RV survey for close companions to very 
low-mass objects in the Chal star-forming region based on 
UVES/VLT spectra. The survey is sensitive to companions at 
orbital distances of 3 AU and smaller (averaged over all objects) 
with a detection rate of 50% or more, as shown by a Monte- 
Carlo simulation. For BD/VLM binaries with a mass ratio close 



to unity {q >0.8), the survey is sensitive to even larger sep- 
arations. This is a significant extension to larger orbital sepa- 
rations compared to previous results of this survey (Joergens 
2006a), as well as compared to other RV surveys of young 
BD/VLMS (Kurosawa et al 2006; Prato 2007a; Maxted et al. 
2008), which cover separations <0.3 AU, and to other surveys of 
field BD/VLMS (Reid et al 2002; Guenther & Wuchterl 2003; 
Basri & Reiners 2006), which cover separations <0.6 AU (com- 
parison based on q between 0.2 and 1, 3.3 o" detection, 50% de- 
tection probability). Only the survey by Basri & Reiners (2006) 
has a comparable time base, however, the survey presented here 
is, nevertheless, superior in terms of detection efficiency mainly 
because of smaller RV errors. Thus, for the first time, the bi- 
nary frequency of BD/VLMS is probed for the whole separation 
range <3 AU with a sensitivity to low-mass ratio systems. 

Two objects were identified as spectroscopic binaries based 
on significant RV variability. This corresponds to an observed 
binary fraction of 18^j2 % for the whole sample of eleven BDs 
and (very) low-mass stars (M< 0.25 Mq) and lO^g^ % for the 
subsample of ten BD/VLMS (M< 0.1 Mq), respectively. The de- 
tected binaries are (i) the BD/VLMS ChaHcK8 (M5.75-M6.5), 
which has a substellar RV companion in a ~1 AU orbit, as pre- 
viously discovered (Joergens 2006; Joergens & Miiller 2007). 
Cha Hcc 8 is likely the lowest mass RV companion found so far 
around a BD/VLMS. (ii) The low-mass star CHXR 74 (M4.25- 
M4.5), which has a spectroscopic companion in a relatively 
long -period orbit (presumably >12 years), as found by previous 
observations (Joergens 2006) and confirmed here based on new 
RV data. With such a relatively long period, CHXR 74 is com- 
ing into reach for being one of the first very young low-mass 
spectroscopic binaries that can be spatially resolved. 

These spectroscopic binaries appear to have orbital periods 
of at least a few years, i.e. orbital separations of 1 AU or greater: 
(i) an RV orbit solution shows that Cha Ha 8 has a companion in 
an orbit of more than four years (Joergens & Miiller 2007); (ii) 
the companion of CHXR 74 has very likely an orbit of more than 
twelve years (Sect.|5]). Thus, while the rate of BD/VLM binaries 
at <3 AU is found to be 10% (18% including CHXR 74), at sep- 
arations <1 AU, it might be less than 10% (0/11 for the whole 
sample and 0/10 for the BD/VLMS subsample). There were no 
signs of the presence of shorter period companions around the 
targets. This is noteworthy because those cause stronger RV 
signals making them easier to detect. In fact, two out of the 
four BD/VLM spectroscopic binaries with determined orbital 
parameters (including Cha Ha 8) have periods of only a few days 
(Basri & Martin 1999; Stassun et al. 2006). 

There is no evidence in the cross correlation peaks of a 
double-lined nature of the detected binaries. This hints at a low 
rate of double-lined spectroscopic binaries in the sample. 

In this survey, there are no large velocity amplitudes 
recorded (Fig. [13]), as would have been expected for (nearly) 
equal mass binaries (e.g. a few tenths of kms ' for the BD 
binary PPl 15, Basri & Martin 1999). The detected significant 
RV variability occurs only on small scales: Cha Ha 8 has an 
RV semi-amplitude of 1.6kms ' (Joergens & Miiller 2007) and 
a half peak-to-peak RV difference for CHXR 74 is 2.4kms 
Note that an only moderate RV precision of e.g. 1.4kms"\ as 
in the survey by Basri & Reiners (2006), would have prevented 
the detection of these binaries. The results hint at the possibility 
that there exist more low-mass ratio systems among BD/VLMS, 
at least at <3 AU, than anticipated from an extrapolation from 
the mass-ratio distribution derived based on direct imaging sur- 
veys, which peaks at ^ ~ 1 and declines rapidly for low q val- 
ues (e.g. Burgasser et al. 2007). On the other hand, the small 
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RV amplitudes can also be attributed to observations at low or- 
bital inclinations. It is, however, interesting that a preference for 
small RV amplitudes is also seen in RV surveys of BD/VLMS 
in other star-forming regions, namely in USco and p Oph (0.2- 
0.5 km s"', Kurosawa et al. 2006, cf. also footnote 5) and in 
Taurus, Ophiuchus, and TW Hya (1.5-2.0 km s ', Prato 2007a). 

Current spectroscopic surveys of very low-mass objects dif- 
fer widely in the studied parameter ranges, as orbital separa- 
tions, (primary and) companion masses, and environments/ages. 
Within the (sometimes large) statistical errors, the determined 
binary frequencies are largely consistent and range between 
<7.5% and 24%: 10!^*^ % for ten BDA'LMS in Chal at <3 AU 
(this work), 24+1^ % for 17 BD/VLMS in USco and pOph at 
<0.3AU (Kurosawa et al. 2006fl 17% for 18 BD/VLMS in 
Taurus, Ophiuchus, and TW Hya at <0.3AU (Prato 2007a), 
<7.5% for 51 BD/VLM member candidates of cr and A Orionis 
at <0.3 AU (Maxted et al. 2008), 12+'' % for 25 field and cluster 
dwarfs at <0.6 AU (Guenther & Wuchterl 2003), and 11+^ % for 
53 (very) low-mass field stars and BDs at <0.2 AU (<l-2 AU for 
q >0.8) (Basri & Reiners 2006). 

By combining my work in Chal with other RV surveys of 
young regions, an overall BD/VLM binary fraction at an age 
of a few million years can be derived. With the exception of 
Chal, all RV surveys of young BD/VLMS probe only a nar- 
row separation range (<0.3 AU). No companions are found or- 
biting BD/VLMS in Chal at these separations. Thus, the overall 
observed binary fraction at separations <0.3 AU of very young 
BD/VLMS in Chal (0/10), USco and p Oph (4/17, Kurosawa et 
al. 2006), Taurus, Ophiuchus, and TW Hya (3/18, Prato 2007a), 
and o- and A Orionis (0/52, Maxted et al. 2008) is Tl% (7/97). 

Comparing the binary frequency of BD/VLMS in Chal for 
separations <3AU (10;^g^%) with what is determined based 
on direct (adaptive optics and HST) imaging surveys at larger 
separations (>20AU) in the same star-forming region (9^^^%, 
Neuhauser et al 2002; Neuhauser, Guenther & Brandner 2003; 
11+^%, Ahmic et al 2007), one finds consistent values. 

7. Conclusions 

For the first time, the binary frequency of BD/VLMS is probed in 
the orbital separation range 1-3 AU. The separation distribution 
of BD/VLM binaries detected by direct imaging surveys (e.g. 
Burgasseret al. 2007) has a peak around 3-10 AU, which is close 
to the incompleteness limit of these surveys (<2-3 AU for the 
field and <10-15AU for star-forming regions, respectively). It 
seemed, therefore, possible that the largest fraction of BD/VLM 
binaries has separations in the range of about 1-3 AU and still re- 
mained undetected. The BD/VLM binary frequency of 10;^^** % 
determined here at separations of 3 AU and smaller shows that 
this is not the case because the rate of spectroscopic binaries at 
<3 AU does not significantly exceed the rate of resolved bina- 
ries at larger separations (10-30%, e.g. Burgasser et al. 2007). 
Thus, the separation range <3 AU missed by direct imaging sur- 
veys is not adding a significant fraction to the BD/VLM binary 
frequency, and the overall binary frequency of BD/VLMS is in 

The four objects classified as spectroscopic binaries in the survey 
by Kurosawa et al. (2006) display very small RV differences (0.2- 
0.5 km s"'). It has to be further investigated whether it can be excluded 
that they are caused by systematic errors. If interpreted as RV semi- 
amplitudes caused by companions at <0. 1 AU, these companions would 
have masses M2 sin i in the planetary regime (1-2 Mj^p) (see Sect[3]l. In 
both cases, systematic errors or planets, we would not count them as 
binaries yielding a binary fraction of this survey of < 6% (0/17). 



the range 10-30%. The finding of no signs (0/10, i.e. <10%) for 
companions at <1 AU is consistent with a decline in the sepa- 
ration distribution of BD/VLM binaries towards smaller sepa- 
rations starting between 1 and 3 AU. Furthermore, the previous 
finding that the observed mass-dependent decrease in the stellar 
binary frequency (e.g. Delfosse et al. 2004; Shatsky & Tokovinin 
2002; Kouwenhoven et al. 2005) continuously extends into the 
BD regime (e.g. Burgasser et al. 2007) is confirmed here for sep- 
arations <3 AU. This is consistent with a continuous formation 
mechanism from low-mass stars to BDs. 

These results are based on a relatively small sample. By com- 
bining this work in Cha I with surveys of other young regions, 
which are restricted to separations <0.3 AU, an overall observed 
binary fraction of very young BD/VLMS can be determined with 
better statistics but with a more limited separation range. In this 
way, a binary frequency of 7'^^% (7/97) is found at <0.3 AU. In 
the near future, the sample of the survey in Chal will be sub- 
stantially enlarged, allowing the re-investigation of the binary 
frequency at <3 AU of BD/VLMS in Cha I with significantly im- 
proved statistics for a uniformly obtained data set. Current and 
future observational efforts are, in addition, directed towards de- 
termination of orbital parameters of the detected binary systems, 
including RV follow-up and high-resolution imaging/astrometry. 



Table 1. 
Chal. 
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Ohiert 


Date 


HJD 


RV 

[kms"'] 


[km s"'] 


CliaHal 


2006 Jan 16 


2453751.74004 


16.206 


0.39* 


Cha Hff 2 


2006 Jan 16 


2453751 79554 


16.319 


0.25 * 


diaHaS 


2006 Jan 15 


2453750.81659 


16.339 


0.70* 


CliaHa5 


2006 Jan 15 


2453750.84774 


15.697 


0.03* 


CliaHae 


2006 Jan 16 


2453751.82994 


16.630 


0.07* 


CliaHa7 


2006 Jan 15 


wrong pointing 






CliaHQ'12 


2006 Jan 17 


2453752.70145 


15.831 


0.05* 


B34 


2005 Nov 25 


2453699.83976 


15.687 


0.10* 


CHXR74 


2005 Mar 21 

2006 Apr 10 
2006 Jun 15 


2453450.58923 
2453835.63456 
2453901.52267 


18.185 
17.982 
19.063 


0.06* 

0.10 

0.10 



" Error of the relative RVs. In most cases (marked with an asterisk) 
this is the sample standard deviation derived from two consecutive mea- 
surements. 



Table 2. Characteristics of observations and;i'^ probabilities. 



Object 


ARV^ 
[kms~] 


RV 

[kms-'] 


no. obs_^ 


Min. At 

[days] 


Max. At 

[days] 


P 


Cha Hal 


0.24 


16.27 ± 0.52 


3 


20 


2113 


8.3 xlO-' 


Cha Ha 2 


0.15 


16.25 ± 0.49 


3 


20 


2113 


8.7 xlO-' 


Cha Ha 3 


0.99 


14.86 ±0.89 


3 


19 


2111 


5.6 xlO-- 


Cha Ha 4 


0.21 


14.85 ± 0.43 


11 


1 


701 


1.0 


Cha Ha 5 


0.13 


15.59 ±0.48 


3 


19 


2111 


8.7 xlO-' 


Cha Ha 6 


0.28 


16.52 ±0.56 


3 


19 


2112 


6.3 xlO-' 


Cha Ha 7 


0.58 


17.09 ± 0.98 


2 


19 


19 


1.5 xlO"' 


Cha Ha 8 


1.38 


16.62 ± 0.67 


11 


3 


2542 


0.0 vari. 


Cha Ha 12 


0.96 


15.29 ±0.93 


3 


20 


2113 


5.6 xlO-^' 


B34 


0.07 


15.74 ± 0.42 


6 


6 


2083 


1.0 


CHXR74 


2.39 


16.71 ±0.81 


13 


1 


2285 


0.0 vari. 



" half peak-to-peak RV differences 

* the error of the weighted mean takes an uncertainty of 400 ms"' 
into account for the zero point of the velocity 

includes data from Table[T] Joergens (2006a), and Joergens & 
Muller (2007) 
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Fig. 3. RV measurements of Cha Ha 7 and Cha Hor 12 (top to 
bottom). 
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Fig. 1. RV measurements of Cha Ha 1, Cha Ha 2, Cha Ha 3 

(top to bottom). 
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Fig. 4. RV measurements of B34 and CHXR74 (top to bot- 
tom). 
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Fig. 2. RV measurements of Cha Ha 4, Cha Ha 5, Cha Ha 6 

(top to bottom). 



Joergens: Binary frequency of very young brown dwarfs at separations <3 AU 



9 



10 ' 



> 5" 



3 

-log p 



Fig. 5. Histogram of the probability p for fitting the ob- 
served (relative) RV values of the BD/(V)LMS in Cha I studied 
here with a constant function. Objects with p < 10"^ are shown 
in the right-most bin, ChaHorS and CHXR74. Based on the 
threshold p < 10"^, they are classified as binaries. Overplotted is 
the expected distribution for nonvariable objects given the esti- 
mated uncertainties. It shows that the error estimates are reliable. 
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Fig. 6. Detection probability as function of semi-major axis. Based on a Monte Carlo simulation of RV measurements for ran- 
domly selected binaries using the measurement errors and time separations of the real observations. From top to bottom: Cha Ha 7; 
average of ChaHorl and ChaHar3; average of ChaHa2, ChaHa4, ChaHorS, ChaHar6, ChaHaS, and ChaHal2; B34; and 
CHXR74. As indicated, some curves are the average detection probability for two or more objects because observations have 
similar quahty for those. 
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Fig. 7. Dependence of detection probability on assumption about inclination. Both curves are calculated for Mi =0.07 M©, 

crxv=0.29 km s ^ fofc5=|0,20d,2000dj, and (7 selected randomly from |0.2,1]. Assuming a fixed mean inclination value (dotted line), 
as in Basri & Reiners (2006), slightly overestimates the detection efficiency compared to the realistic case of random orientation 
(soUd fine). For the fixed inclination case, the mean value of the inclination distribution of 57.3° is used. 



12 Joergens: Binary frequency of very young brown dwarfs at separations <3 AU 




Fig. 8. Dependence of detection probability on number of observations. Both curves are calculated for Mi =0.07 M©, 
o'i?v=0.29kms~', a maximum time span of 2000 d, q selected randomly from [0.2,1], and random orientation. An observing sched- 
ule based on three observations with foi,5=[0,20d,2000d] (solid line) provides a superior detection efficiency compared to the case of 
only two observations separated by 2000 d (/oij=[0,2000d], dotted hue). Thus, third-epoch observations between the observations 
with the largest time difference can have a significant effect on the probability distribution. 
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Fig. 9. Comparison of detection probability of different BD/VLMS RV surveys. All curves were calculated for q selected rati- 
domly from [0.2,1], random orientation, variability threshold of p < 10"^^. Solid (red) line: this work (calculated for Mi-O.Qil M^, 
o-Rv=0.29kms-i, foi,=[0,20d,2000d]); dotted (black) line: Basri & Reiners 2006 (Mi=0.15 M^, o-rv=1.3 kms"', fofe,,=[0,1650d]); 
dashed (blue) Hne: Guenther & Wuchterl 2003 (Mi=0.1Mo, o-Rv=0.64kms"', foi,=[0,31 d,85 d]); dash-dotted (green) Hne: 
Kurosawa et al. 2006 (Mi =0.07 Mq, crRV'=0.42kms"', foi(,j=[0,20d]). See online version for a color figure. 
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Fig. 10. Same as Fig.|9]but with logarithmic y-axis to facilitate comparison with other publications. See online version for a 
color figure. 
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Fig. 11. Comparison of detection probability of different BD/VLMS RV surveys for high-mass ratio binaries. Same as Fig.|9] 
but for q selected randomly from [0.8,1]. See online version for a color figure. 
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Fig. 12. Same as Fig.[TT]but with logarithmic y-axis to facilitate comparison with other publications. See online version for a 
color figure. 
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Fig. 13. RV differences vs. object mass. Plotted are half peak- 
to-peak differences in the observed RVs. Each data point is 
labeled with the corresponding object name, the numbers de- 
note the ChaHo' objects. The RV variable objects (CHXR74, 
ChaHo-S) have ARV > 1.4km s ', while the majority of RV 
constant objects have ARV <0.3-0.6kms"', with the exception 
of Cha Ha 3 and Cha Hot 1 2, which have ARV ~ 1 km s ' . 
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Fig. 14. RV differences for BD/VLMS in USco and pOpli 

based on data from Kurosawa et al. (2006). ARV as in Fig.[T3] 
Objects classified as RV variables by these authors are denoted 
with their names. The two data points with the largest RV differ- 
ence correspond to GY310 and UScolOl, which are classified as 
RV constant. 
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